Interface Characteristics at an Organic/Metal Junction: Pentacene on Cu Stepped Surfaces

The adsorption of pentacene on Cu(221), Cu(511) and Cu(911) is investigated using density functional theory (DFT) with the self-consistent inclusion of van der Waals (vdW) interactions. Cu(211) is a vicinal of Cu(111) while Cu(511) and (911) are vicinals of Cu(100). For all the three surfaces, we found pentacene to prefer to adsorb parallel to the surface and near the steps. The addition of vdW interactions resulted in an enhancement in adsorption energies, with reference to the PBE functional, of around 2 eV. With vdWs inclusion, the adsorption energies were found to be 2.98 eV, 3.20 eV and 3.49 eV for Cu(211), Cu(511) and Cu(911) respectively. These values reflect that pentacene adsorbs stronger on (100) terraces with a preference for larger terraces. The molecule tilts upon adsorption with a small tilt angle on the (100) vicinals (about a few degrees) as compared to a large one on Cu(221) where the tilt angle is found to be about 20o. We find that the adsorption results in a net charge transfer to the molecule of ~1 electron, for all surfaces.Comment: 11 pages, 4 figure

Electronic Structure of the c(2x2)O/Cu(001) System

The locally self-consistent real space multiple scattering technique has been applied to calculate the electronic structure and chemical binding for the c(2x2)O/Cu(001) system, as a function of $d_{O-Cu1}$ -- the height of oxygen above the fourfold hollow sites. It is found that the chemical binding between oxygen and copper has a mixed ionic-covalent character for all plausible values of $d_{O-Cu1}$. Furthermore, the electron charge transfer from Cu to O depends strongly on $d_{O-Cu1}$ and is traced to the variation of the long-range electrostatic part of the potential. A competition between the hybridization of Cu1-$d_{xz}$ with O-$p_x,p_y$ and Cu1-$d_{x^2-y^2}$ with O-$p_z$ states controls modification of the electronic structure when oxygen atoms approach the Cu(001) surface. The anisotropy of the oxygen valence electron charge density is found to be strongly and non-monotonically dependent on $d_{O-Cu1}$.Comment: 14 pages, 7 figures, 1 tabl

The crossover from collective motion to periphery diffusion for 2D adatom-islands on Cu(111)

The diffusion of two dimensional adatom islands (up to 100 atoms) on Cu(111) has been studied, using the self-learning Kinetic Monte Carlo (SLKMC) method [1]. A variety of multiple- and single-atom processes are revealed in the simulations, and the size dependence of the diffusion coefficients and effective diffusion barriers are calculated for each. From the tabulated frequencies of events found in the simulation, we show a crossover from diffusion due to the collective motion of the island to a regime in which the island diffuses through periphery-dominated mass transport. This crossover occurs for island sizes between 13 and 19 atoms. For islands containing 19 to 100 atoms the scaling exponent is 1.5, which is in good agreement with previous work. The diffusion of islands containing 2 to 13 atoms can be explained primarily on the basis of a linear increase of the barrier for the collective motion with the size of the island

Island Size Selectivity and island-shape analysis during 2D Island Coarsening of Ag/Ag (111) Surface

In our earlier study of Ag island coarsening on Ag(111) surface using kinetic Monte Carlo (KMC) simulations we found that during early stages coarsening proceeds as a sequence of selected island sizes resulting in peaks and valleys in the island-size distribution and that this selectivity is independent of initial conditions and dictated instead by the relative energetics of edge-atom diffusion and detachment/attachment processes and by the large activation barrier for kink detachment. In this paper we present a detailed analysis of the shapes of various island sizes observed during these KMC simulations and show that selectivity is due to the formation of kinetically stable island shapes which survive longer than non-selected sizes, which decay into nearby selected sizes. The stable shapes have a closed-shell structure - one in which every atom on the periphery having at least three nearest neighbors. Our KMC simulations were carried out using a very large database of processes identified by each atom's unique local environment, the activation barriers of which were calculated using semi-empirical interaction potentials based on the embedded-atom method.Comment: 17 pages, 11 figure

SLKMC-II study of self-diffusion of small Ni clusters on Ni (111) surface

We studied self-diffusion of small 2D Ni islands (consisting of up to 10 atoms) on Ni (111) surface using a self-learning kinetic Monte Carlo (SLKMC-II) method with an improved pattern-recognition scheme that allows inclusion of both fcc and hcp sites in the simulations. In an SLKMC simulation, a database holds information about the local neighborhood of an atom and associated processes that is accumulated on-the-fly as the simulation proceeds. In this study, these diffusion processes were identified using the drag method, and their activation barriers calculated using a semi-empirical interaction potential based on the embedded-atom method. Although a variety of concerted, multi-atom and single-atom processes were automatically revealed in our simulations, we found that these small islands diffuse primarily via concerted diffusion processes. We report diffusion coefficients for each island size at various tepmratures, the effective energy barrier for islands of each size and the processes most responsible for diffusion of islands of various sizes, including concerted and multi-atom processes that are not accessible under SLKMC-I or in short time-scale MD simulations

Atomistic studies of thin film growth

We present here a summary of some recent techniques used for atomistic studies of thin film growth and morphological evolution. Specific attention is given to a new kinetic Monte Carlo technique in which the usage of unique labeling schemes of the environment of the diffusing entity allows the development of a closed data base of 49 single atom diffusion processes for periphery motion. The activation energy barriers and diffusion paths are calculated using reliable manybody interatomic potentials. The application of the technique to the diffusion of 2-dimensional Cu clusters on Cu(111) shows interesting trends in the diffusion rate and in the frequencies of the microscopic mechanisms which are responsible for the motion of the clusters, as a function of cluster size and temperature. The results are compared with those obtained from yet another novel kinetic Monte Carlo technique in which an open data base of the energetics and diffusion paths of microscopic processes is continuously updated as needed. Comparisons are made with experimental data where available

Island Size Selectivity during 2D Ag Island Coarsening on Ag (111)

We report on early stages of submonolayer Ag island coarsening on Ag(111) surface at room temperature ($300$ K) carried out using realistic kinetic Monte Carlo (KMC) simulations. We find that during early stages, coarsening proceeds as a sequence of selected island sizes creating peaks and valleys in the island size distribution. We find that island-size selectivity is due to formation of kinetically stable islands for certain sizes because of adatom detachment/attachment processes and large activation barrier for kink detachment. In addition, we find that the ratio of number of adatom attachment to detachment processes to be independent of parameters of initial configuration and also on the initial shapes of the islands confirming that island-size selectivity is independent of initial conditions.These simulations were carried out using a very large database of processes identified by their local environment and whose activation barriers were calculated using the embedded-atom method

Phonons of Metallic Vicinal Surfaces

We present an analysis of the vibrational dynamics of metal vicinal surfaces using the embedded atom method to describe the interaction potential and both a real space Green's function method and a slab method to calculate the phonons. We report two main general characteristics : a global shift of the surface vibrational density of states resulting from a softening of the force field. The latter is a direct result of the reduction of coordination for the different type of surface atoms; and an appearance of high frequency modes above the bulk band, resulting from a stiffening of the force field near the step atom. The latter is due to a rearrangement of the atomic positions during the relaxation of the surface atoms yielding a large shortening of the nearest neighbor distances near the step atoms.Comment: 6 figures, to appear in Sur. Sci. proceedings of VAS1